ON THE DIFFUSIVE STRESS RELAXATION FOR 
MULTIDIMENSIONAL VISCOELASTICITY 

DONATELLA DONATELLI AND CORRADO LATTANZIO 

Abstract. This paper deals with the rigorous study of the diffusive 
stress relaxation in the multidimensional system arising in the mathe- 
matical modeling of viscoelastic materials. The control of an appropriate 
high order energy shall lead to the proof of that limit in Sobolev space. 
It is shown also as the same result can be obtained in terms of relative 
modulate energies. 



1. Introduction 

The aim of this paper is the study of the diffusive relaxation limit for a 
model in multidimensional viscoelasticity. To this end, we shall consider the 
following semilinear hyperbolic system 

d t F ia - d a Vi = 

d t Vi - d a S ia = (1.1) 
d t S ia - £d a Vi = -\S ia + \T ia (F). 

In (|1.1|) . i, a = 1, 2, 3; F and v are respectively the deformation gradient 
and the velocity, while e > stands for the relaxation parameter. Moreover, 
we assume the stress tensor T to be smooth and, from now on, we sum on 
repeated indices. Putting formally e = in Ijl.lJI . we recover the equilibrium 
relation 

S = T(F)+V x v (1.2) 

and the equilibrium dynamic is hence given by the incomplete parabolic 
system for viscoelasticity 

d t F ia - d a Vi = 
dm - d a T ia (F) = fid a d a v. 

The rigorous proof of this diffusive relaxation limit will be showed in the 
framework of Sobolev spaces, as long as the reduced system (|1.3|) admits 
(regular) solutions (Section |2J). To perform this task, we follow the approach 
already used in J^] to analyze a singular Euler-Poisson approximation for 
the incompressible Navier-Stokes equations and previously in [S] for singular 
perturbations of general first order hyperbolic pseudo-differential equations. 
More precisely, via standard energy estimates, we shall obtain a suitable 
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bound in for the differences v — v, F — F and S — S, for well-prepared 
initial data, namely when they approach equilibrium in that space. This 
result will imply in particular the existence of (classical) solutions of 
in an e-independent time interval [0, T] and their convergence toward the 
solutions of in that interval. Moreover, we shall recast this convergence 
result also in terms of relative modulated energy techniques (Section 12. 

The results contained in this paper are part of a more general project 
which intend to study connections between (semilinear) system of conser- 
vation laws with diffusive source terms and (possible degenerate) parabolic 
systems. The mathematical study of the these connections starts from the 
papers of Kurtz [Sj and McKean [Tl|, where the authors introduced for the 
first time a parabolic scaling for hyperbolic systems, to put into evidence 
their diffusive behavior. Afterwards, this scaling has been extensively used 
in the analysis of hyperbolic-parabolic relaxation limits, both for weak so- 
lutions, by means of compensated compactness techniques (among others, 
see ^3 |7j and the references therein), and for classical solutions (see, for 
instance, |H1 El)- It is worth to observe that, as for the one-dimensional 
model treated in , in the present case the equilibrium system turns out to 
be incompletely parabolic, which forces ourselves to study the relaxation limit 
in the context of regular solutions, because standard compensated compact- 
ness techniques do not apply. In this framework, in addition to the already 
mentioned paper [2] , further results dealing with diffusive relaxation toward 
a degenerate parabolic limit are confined to the case of a class of BGK-type 
approximations to a Cauchy problem for multidimensional degenerate scalar 
parabolic equations p. 

Finally, we point out that, as in 0, the relaxation limit in the present 
case has also a physical interpretation in terms of mathematical models in 
the study of viscoelastic materials More precisely, it can be viewed 

as the passage from the viscosity of the memory type to the viscosity of the 
rate type. Indeed, at the equilibrium, the stress-strain response S is given 
by P-2jl . while in S can be recast as follows 



In this section we study the relaxation limit e | for system p.l|) in 
the framework of Sobolev spaces, showing the convergence of its solutions 
toward the solutions of (|1.3|) . as long as the latter exist. 

Remark 2.1. The local well-posedness in Sobolev spaces for large data of 
the model 1)1.3(1 is guaranteed by the results on general weakly parabolic 
systems jlj Theorem 1 and Remark 2]. To recover the lack of parabolicity 
of the model, the authors proposed in the paper appropriate conditions on 




that is, the viscous effect comes from a memory term. 



2. Energy estimates and trend to equilibrium 
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the first order terms, which must be controlled by the diffusion term. These 
sufficient conditions leading to the local existence result reduce in our case 
to 
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Ra-Aa(F) T )v a \ <C7x^K| 2 , (2.1) 



a=l a=l 



for any F E Mat 3x3 ~ R 9 and for any v a G M 3 , a = 1, 2, 3. In (J2HJ), the 
matrices R a , A a (F) T G Matg x3 are defined by 

'<Wsx3\ (Ti a {Ff 
Ra = I hahxS , A a {F) = V F T 2a (F) 



a(F) Y 



V <5 3Q /3X3/ V r 3a(^)/ V i/3 A 

for a = 1, 2, 3. Hence condition 1)2. 1|) is fulfilled if and only if 

V F T(F) < n, (2.2) 



J/3 



for any F under consideration. Condition 1)2.2(1 has been considered already 
in in the one-dimensional case to prove global existence and convergence 
of relaxation limit in Sobolev norms for large solutions. Moreover, it stands 
for the subcharacteristic condition assumed in to control the relaxation 
limit in the hyperbolic-hyperbolic regime. In the multidimensional case, 
1)2.2(1 gives local existence of smooth solutions to 1)1.3(1 . thanks to while 
global smooth solutions can be constructed for small perturbation-type ini- 
tial data. 

Let (F,v) T be the solution of <JOJ), belonging to C 1 ([0, T]; fF(R 3 )) for 
any j > | + 2 and define S by 1)1.2(1 . Then the differences u = u — u, 
F = F = F and S = S — S verify the following semilinear system 

' d t F ia - d a Vi = 
dm - d a S ia = 

d t S ia - £d aVi = -\S ia + \ (T ia (F + F) — T ia (F)) - d t S ia . 

(2.3) 

We rewrite (|2.3|) in vectorial notations as follows 

W t + A e a d a W = R(W;W), (2.4) 



where 



W = (F, S) T , W = (F, v, S) T , 

/0gx9 — -Ra Og x g\ 

A e a = 3x3 3x3 -C € Mat 2 ix2i, 



Og x g 


—Ra 


Og x g 


3 x3 


3 x3 




Og x g 




0gx9 

















7 5 ia + | (T ia (F + F) — T ia (F)) - d t S it 
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and system IJ2.4JI is coupled with initial data 

W^(x) = (F£(x)-F(x,0),v e o (x)-v(x,0),S e (x) - S(x,0)f. 

In the next lemma, we show system Q2.4|l is symmetrizable and hence it is 
equipped with a positive definite energy, at least for e ^ 1 (see [Sj for the 
same analysis in 1-D). 

Lemma 2.2. For e<l, the matrix 

(^9x9 0gx3 —^9x9 

3 x9 (f " 1) / 3 X3 3x9 
—-^9x9 09x3 -^9x9 

defines a positive definite symmetrizer of j2.4\ ) and 

E € (W) = (B e W,W) L 2 { 
is a positive definite energy for solutions of j2.4[ l- 
Proof. A direct computation shows 

(0gx9 0g X 3 0g X 9 

3 x9 3x3 - (1 ~ 1) Rl 

0gx9 — — l) Ra 0g x g 

that is, B e symmetrizes (j2.4|) . Moreover, 

E?(W) = J (^|F| 2 - 2F a S a + (~ - l) \v\ 2 + \S\ 2 ) dx 

>\ [ (^\F\* + %\* + \Sf)dx, 
2 J R 3 \e e J 

for e < 1, say e < \. □ 

The above result implies the energy eE e (W) is equivalent to the (square 
of the) 1? norm of (F(-,t),v(-,t), ^eS(-,t)) T , for any t > and for e < 1. 
Thus we define the following high-order energy: 

£\t) = eE^W), (2.5) 

W<3 

where 7 is a multi-index. Hence we shall establish our convergence result in 
terms of the norm defined in (|2.5|) , which is equivalent for e sufficiently small 
to the (square of the) H 3 norm of the vector (F(-, t), v(-, t), y/eS(-, t)) T . 

Theorem 2.3. Let (F, v) T be a solution of fl.fy) belonging to the space 
C 1 ([0, T]; ff 4 (R 3 )) and define S by ll.ty) . Let us assume the initial data 
(F$,v € ,S^) T E H 3 (R 3 ) for E2J) verify 

5 £ (0) = e ^(W) - 0, 

ItI<3 
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as e I 0. Then there exist eo > and a constant Ct > such that for 
any < e < eo and any t € [0, T\, there exist a (unique strong) solution 
(F(;t),v(-,t),V~eS(-,t)) e F 3 (M 3 ) of [Hi) verifying 

S e (t) = eE'id^W') < C T (e 2 + £ e (0)). (2.6) 

Proof. Prom classical results of local well-posedness of Sobolev solutions for 
symmetrizable hyperbolic systems (see, for instance, |13l 1181 IT7] ). we know 
there exists a maximal time T e > such that, for any t G [0, T e ), 

£ e (t) < M\ (2.7) 

where the value M e , decaying to zero as e j 0, will be chosen later. The result 
of the theorem is then equivalent to prove T e > T, which shall be obtained 
by showing equality in (|2.7|) cannot be achieved for T e < T, provided M e is 
properly chosen jT^]. 

To perform the energy estimate needed to prove (|2.6|) . we apply the sym- 
metrizer B e to the partial derivative of the nonhomogeneous term R(W\ W), 
and we obtain 

1 \d^S ia - ±d 7 (T ia (F + F) - T ia (F)^ + d t d 7 S ia ^ 
B e d y R{W;W) = 

K -^S ia + ±c> 7 (T ia (F + F) — T ia (F)) - d t d 7 s iaj 

We start by computing the time derivative of the energy E e {d 1 W). By 
integrating by parts and by taking into account Lemma 12.21 we get 

^E'id^W) = 2{B e d 1 R{W-W) 1 d~ j W) L 2 { 



2 

+ (r ia {F + F) - T ia {F)^j , d 1 S ia - d 7 F ia ) i2 

2 ~ 
||9 7 5||£,2( R 3) + (dtdrySiajd^Fia — d^Si a ) L 2 



€ 

= h+h + h + h (2.8) 
Now we estimate separately each one of the terms I\, I2, In- For I\ we have 



2 12, 

h = -(9-ySi a ,d 7 Fia)L2(M3) < ^-||^7'S'|| i 2( R 3) + -||9 7 F|| . (2.9) 

The term I2 can be controlled as follows 

h = (T ia (F + F) - T ia (F^j , d 7 S ia - d 7 F 4Q ) L 2 

= -{d 1 (v k TI Ti a (F)F] CV \ ,djSi a — d^Fi a ) L 2 



€ 

+ _ (^7 [Tia(F + F) — T ia (F) — Vkr 1 Ti ol (F)F kr] j ,d 1 S ia — d 1 F ia ) L 2 
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- ^H^t^Hl^r 3 ) + ~II^7-^IIl 2 (r 3 ) + ~H^r \^k v Ti a {F)F kri ) \\ L2 
+ ^||d 7 n d$J duVlxTiaiuW + ^/b^c) lli 2 (R3), (2-10) 
where we used the notations 

V T (F) 9Tla{F) V 2 T (F) - 92Tla{F) 

Finally we estimate J4, 

I4 = 2(dtd y Sia, d 1 Fi a — d 1 Si a ) L 2( 

< 2e||<9 t a r S'|| 2 : 2( K 3) + -[|9 7 5[|^2( R 3) + -||9 7 -F[|i2( K 3). ( 2 -ll) 
By using (|23j) . (t2~TUl) and (l2~TTT) in (JUJ) we get 

= 2e||a t a 7 5||| 2(R3) - ^||5 7 5||i 2(M3) + ^ 7 F||i 2(R3) 
+ ^||d 7 (£d0^dujVl vlc T ia (uOF + F)6F kv F lc } || 2 2 

+ ~ll^7 {y kr)Tiot{F)F kr ^j |||2( R 3) (2.12) 

Now we multiply (|2.12|) by e we sum in 7, I7I < 3 and by taking into account 
Sobolev inequalities we end up with 

±F(t) + ]\\Sf H s (m < c(\\F\\ L ^(t) + c(\\F\\l~, \\F\\ L ~)(£'(t) 2 



dt 



+ £%ty) + 2^\\d t S\\U m . (2.13) 



To simplify notations, we denote with ~ all constants depending on W and its 
derivatives. Fixing M t < 1, in particular we have Hi^Hioo^s) < 1 and hence 

there exists a constant C2 > such that c(||F[|£,oo) + c(||F||/,oo , [|F||_l=x>) < cjj. 
Then, from (|2.13j) it follows 

— £ e (t) < cie 2 + c 2 £ e (t), for any t G [0, T e ] 

and the Gronwall Lemma implies 

£ e {t) < (£ e (0) + cie 2 t) e^ 2t for any t G [0,T e ]. (2.14) 

If we choose M e = (<? £ (0) + cie 2 T) /z , we see that, for e sufficiently small, 
we cannot reach inequality in (|2.7j) for T e < T. This proves that T e > T 
and thus that (|2.14|) is valid on [0, T], which conclude the proof. □ 
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2.1. Relative modulated energy approach. The control of the relax- 
ation limit contained in the above theorem can be recast in terms of relative 
modulated energy techniques, already used in |2j for semilinear relaxation 
approximation of incompressible Navier-Stokes equations and in |19| II 1 j in 
the hyperbolic-hyperbolic stress relaxation for elasticity with memory. To 
illustrate this method, we shall obtain here the L 2 control of the difference 
between the equilibrium (F,v), solution of (jOj) . and its relaxation approx- 
imation (F,v). To this end, let us observe that we can eliminate in system 
(|l.lj) the off-equilibrium variable S to obtain 

1 t F ia - d a Vi = ^ 
\d t Vi - d a T ia (F) = fid a d a Vi - edfoi. 

Then, let (F, v) be a smooth, uniformly bounded solution of (|2.15j) for any 
t 6 [0, T] and let us denote with T the corresponding bound on the charac- 
teristic speed of the relaxation approximation, that is 

VfT(f) < n, 

for any F under consideration. For A > 1 arbitrary, we define our relative 
modulated energy as follows: 

W rm = i (\v - v\ 2 + \F- F| 2 ) + e(vi - Vi)d t (vi - + -e 2 X\d t (v - v)\ 2 

+ -e\n\V a (y -v)\ 2 + e\d a {vi - Vi)(T ia (F) - T ia {F)), 
with associated flux given by 

Qa,rm = {v% ~ Vi){T ia {F) - T ia (F)) + flfa - Vi)d a (Vi - Vi) 

+ e\fj,d t (vi - Vi)d a (vi - Vi) + e\d t (vi - Vi)(T ia {F) - T ia {F)). 

Such an energy is obtained modulating the standard energy estimates of 
(jl.3f) by higher order contributions of acoustic waves, to take advantage of 
the dissipation coming from the relaxation term. The above energy verifies 

d{K rm ~ d a Q a , rm + (fJ-\V a (v - v)\ 2 - e\d a (Vi - Vi)V jpTiatFjdpiVj - Vj)) 

+ e(X-l)\d t (v-v)\ 2 

= d a {vi - Vi){F ia - F ia - (T ia {F) - T ia {F))) 
- ed 2 Vi(vi - Vi) - e 2 \d 2 Vid t (vi - Vi) 
+ e\d a {v t - Vi) (y jP T ia (F) - V^T ia {F)) d t F jP 

:=K e . (2.16) 

Before proving (|2.16() . let us emphasize that, for e <C 1, say e < min{^, 1}, 
and for A > 1 properly chosen, 

(m|V«(« - v)\ 2 - eXd a (vi - Vi)V jf} T ia {F)dp{vj - %)) > Ci|V a (tj - v)\ 2 , 
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C 2 <p e (t) > I Hrmdx > C 3 <p*(t), 

where 



<p'(t) := [ (\v-d\ 2 + \F-F\ 2 + e 2 \d t (v-v)\ 2 + e\V a (v-v)\ 2 )dx 

and Ci, C2, C3 > depend on fj, and T and not on e. Moreover, using again 
Schwartz inequality, we control the error terms in H2.16j) to obtain 

\n e \ < (e 2 ^i + K 2 n rm ) + yI v «( € " ^l 2 ' 

with i^i, > depending only on fi, T and the equilibrium solution (F,v) 
and its derivatives and not on e. Hence, we integrate ()2.16j) with respect to 
x and we use Gronwall Lemma to get for any t 6 [0, T] 

^(t)<C T (e 2 + ^(0)), 

that is, the same kind of control of the relaxation limit obtained in Theorem 
I2.3| for well-prepared initial data, that is for <p e (0) — > as e [ 0. 

To prove (j2.16j) . we first observe that the differences F — F and v — v 
verify 

\d t (F ia - F ia ) = d a (vi - Vi) 

\d t (vi - Vi) = d a (T ia (F) - T ia (F)) + fid a d a (vi - v t ) - edf - - edfvi. 

(2.17) 

Then we multiply (|2.17jh by Fi a — Fi a and (|2.17jl o by Vi — Vi and we sum 
over all indices to obtain 

9t - v\ 2 + \F - F\ 2 ) + efa - Vi)d t (vi - 

- d a ((vi - Vi)(T ia (F) - T ia (F)) + ufa - Vi)d a (vi - Vi)j 
+ fi\V a (v-v)\ 2 -e\d t (v-v)\ 2 

= d a ( Vi - Vi)(F ia - Fi a ~ (Ti a (F) - T ia {F))) - ed 2 v t ( Vi - 

(2.18) 

To take advantage of the dissipation coming from the relaxation, we must 
modulate the above relation with an higher order energy, coming from the 
damped wave equation in l|2.17|h . To this end, we fix A > 1 and we multiply 
this relation by e\dt(vi — Vi) and sum over i to obtain 

dt {^e 2 \\d t {v -v)\ 2 + \e\n\V a (v - v)\ 2 

- d a {e\^d t {vi - Vi)d a (vi - Vi)j 

+ e 2 \\d t (v -v)\ 2 - eXdtivi - Vi)d a (T ia (F) - T ia {F)) 
= -e 2 \d t v % d t {vi-Vi). (2.19) 
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We interchange x and t derivatives in the last term in the left of (j2.19|) as 
follows 

- d t (vi - Vi)d a (T ia (F) - T ia (F)) = -d a {vi - v l )d t (T ia (F) - T ia (F)) 

+ d t (d a (vi - Vi){T ia {F) - T ia (F))) ~ d a (fltfa - Vi){T ia {F) - T te (F))) 
= -d a {vi -Vi)S7jpT(F)dp(vj -%) 

- d a (vi - vi) (v j0 T ia (F) - V i/3 T; Q (F)) d t F jf3 

+ d t (d a (vi - Vi){T ia {F) - T ia {F)f) - d a (dt(vi - Vi){T ia {F) - T ia (F))) . 
Finally, using this identity, (|2l9|) and lf2?18|) give 
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